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Macroscopic thermodynamics of equilibrium is constructed for systems obeying 
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^ power-law canonical distributions. With this, the connection between macroscopic 

O 

thermodynamics and microscopic statistical thermodynamics is generalized. This is 
> 

rN complementary to the Gibbs theorem for the celebrated exponential canonical 

distributions of systems in contact with a heat bath. Thereby, a thermodynamic basis is 
provided for power-law phenonema ubiquitous in nature. 
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Statistical mechanics builds an essential bridge between the laws of nature governing 
microscopic dynamics of constituents of matter and its macroscopic behavior. One and 
only example of such a description known to date is the theory of Boltzmann and Gibbs 
characterized by the exponential distributions. Since many systems fall under the sway 
of this theory, there has been no attempt to find alternate possibilities. However, 
nowadays it is widely recognized that many phenomena in nature obey different kinds 
of distributions, i.e., the power-law distributions. We may cite here a few such 
examples: anomalous diffusion [1], vibrating powders [2], fully developed turbulence 
[3], and line shape cumulants in glasses [4]. We point out that even though these are 
typically nonequilibrium phenonema the structures of their distributions persist for 
inordinately long times so that it is natural to consider them as in quasi thermal 
equilibrium. In traditional Boltzmann-Gibbs theory, the canonical distribution function 
is derived from microcanonical ensemble with the principle of equal a priori probability 
(microscopic statistical concept) and is connected with macroscopic thermodynamic 
description by considering the system in contact with a macroscopic heat bath. That it is 
of crucial importance to establish the interrelationship between the two procedures was 
stressed by Tisza and Quay [5] and also by Kubo et al. [6]. 

Recently, it has been shown [7,8] that the power-law distribution can also be deduced 
based on microcanonical ensemble theory, thus providing a foundation to microscopic 
statistical thermodynamics of equilibrium in such a context. It has also been shown 
from the perspective of the (Levy-Gnedenko) generalized central limit theorem [9] that 
this result is a natural extention of corresponding Boltzmann-Gibbs theory obtained 



more than half a century ago by Khinchin [10] using the ordinary central limit theorem. 
However, understanding of the power-law systems in a thermal equilibrium language is 
not complete until a way of realizing their macroscopic thermodynamic description is 
found. 

The purpose of this paper is to present a macroscopic description of systems obeying 
the power-law distributions in equilibrium by reconsideration of contact with a heat 
bath. We establish the interrelationship between microscopic statistical and macroscopic 
thermal equilibria in parallel with the discussion of Boltzmann-Gibbs theory made in 
Ref. [5,6]. Thus, the Gibbsian measure is shown to be appropriately modified here. In 
this manner, a link between the laws on a microscopic scale and the macroscopic 
behavior of matter is extended to a much wider class of systems in nature. 

We begin our discussion by considering a partition of energy between two systems, I 
and II, in thermal contact. These systems have energies E^ and E^^, and the total 
energy is fixed as E - E^ + E^^ . If the state densities of I and II are denoted by ^j and 
^jj, respectively and that of the total system I-i-II by Q. , then we have 

E 

^{E)AE = AE j^^{E^) ^^^{E - E^)dE^ , (1) 



where AE is the width of the shell of constant energy in phase space of the total system. 
The probability of finding the system I in the range of the energy [E^,E^ + dE^) is 
given by [6] 



, , a,{E,) ^JE - E,)AE 
piE, ) dE, = '^ '' "\ '-^ dE, . (2) 



From eq. (1), the normalization of this probability follows obviously. We are interested 
in the most probable partition of energy in conformity with the principle of equal a 
priori probability (microcanonical measure). This is determined by maximizing the 
quantity 

a^{E^)a^^{E - E^)AE dE,. (3) 

In the traditional discussions of identifying the temperature of the system I and deriving 
its probability distribution, one analyzes the maximization procedure by the 
corresponding condition on the logarithm of the quantity in eq. (3). The underlying 
assumption hidden behind this treatment is extensivity (additivity) of entropy. To 
accomodate the power-law distributions of interest in the present work, we relax this 
assumption. From among many possible choices of accomplishing this, here we adopt, 
instead of the ordinary logarithmic operation mentioned above, the "^-logarithmic 
function" defined by [11] 



jc^"^ - 1 

Kix) = -^ -, (4) 

l-q 



where q is a positive constant. This function has the following nonadditive property: 



ln^(xy) = ln^(^) + ln^(3;) + (l-^)ln^(^)ln^(3;). (5) 



The inverse function associated with eq. (4) is the "^-exponential function" 



A^)- 



[l + {l-q)x]"^'"'^ (l + (l-^)x>0) 
(l + (l-^);c<0)' 



(6) 



that is, e^{ln^{x))-x-\n(e^{x)] for positive x. In the limit <? ^ 1, In^(x) and 
e^ix) converge to the ordinary \n{x) and e", respectively. An important point is that, 
just like ln(.x:), ln^{x) is monotonically increasing concave function of x for all 
positive values of q. (It is noted that this particular choice of the pair, ln^(.)c) and e^{x), 
is one of many other possibilities.) Using the ^-logarithmic operation, the maximization 
condition is now written as follows: 

\nXa,{E,)a,,{E-E,)) = max. (7) 

Vanishing of the first derivative with respect to E^ leads to 

1 ^in,(a,(£,)) 



l + {l-q)\(^,{E,)) dE, 



dKi^M)^ (g) 



where equation (5) has been used. If we identify 

5,^ = 5^(a,)^ln^(a,) (/ = I,II) (9) 

as a generalized entropy and define a parameter p. by 



dE, 



then we have 



K = K^P\ (11) 



where we have introduced the notation 

c,^^l + (l-^)5,.^ = a:-^ (12) 

By analogy with the ordinary case {q ^\), the separation constant jS* may be 
interpreted as the temperature of the total system. This can be thought of as the q- 
generalization of the zeroth law of thermodynamics. Similarly, the equilibrium 



conditions associated with the chemical potentials and the pressures can be obtained if 
the particle numbers and system volumes are respectively taken into account in the 
same manner as above. 

If the two systems have different values of q, say q^ and q^^, then the equilibrium 

condition in eq. (11) turns out to be modified to the following form: 



^' ^" =/?*. (13) 



C\q^ '^II^II 



Taking the limit q^^l in this equation, the system I is then regarded as the ordinary 

thermometer measuring the temperature of the system II. Thus, we identify jS* with the 
ordinary inverse temperature ;S, which equates ySy to Cji^ /?. 

To ascertain the above state to be indeed the maximum we are seeking, we need to 
verify that the second derivative is negative. Carrying out such a calculation, we arrive 
at the following expression: 



^' \nX^,[E,)^,{E-E,)) 



dEl 



^(l-^)-^ln[l + (l-^)5,J + ^(l-^)-^ln[l + (l-^)5„,], (14) 



where it is noted that the ordinary logarithm appears on the right-hand side of this 
equation. As in the case of the discussion of ordinary Boltzmann-Gibbs theory, we also 



assume that the second-order derivative of the generalized entropy in eq. (9) is negative. 
Then, we see that each term on the right-hand side of eq. (14) is negative for all positive 
values of q, due to the monotonicity of the logarithmic function. 

Now, let us study the necessary modification of the Gibbs theorem (for the canonical 
measure). For this purpose, we consider the system II to be a heat bath, and therefore 



E,«E,,. 



(15) 



In addition, we take the system I to be in its kih state of energy E^=£j^. The 
probability of finding the system I in such a state is proportional to the number of 
microscopic states of the system II as follows [6]: 



M)' 



a,,{E) 



(16) 



We rewrite this equation in the form 



/(^J 



oc e 



In 



a„(£) 



(17) 



Using the identity 



In, 



'x^ 



yyj 



^[ln,W-ln,(j)], 



(18) 



we obtain 



/(^J-^/7y4^KK(^-£.))-ln,(a,(£))] 



a;-(£) 



(19) 



From eq. (15) with E^ =£4, we may perform the following expansion to the leading 
order of £,: 



f{^.) 



oc e 



1 



a;-(£) 



in^(a,(£))-£, ^y^^ 'U l -in^(a„(£)) 



= ^,(-)S*eJ, 



(20) 



provided that equations (10) and (11) have been used. In the above, the second-order 
term has been neglected because of the assumed large size of the heat bath. Therefore, 
we have 



a,,{E-e,)oca,,{E)eX-P*e,). 



(21) 



This is the modification of the Gibbs measure promised earlier. For q > I, /(f^) in eq. 
(20) for large values of e^, obeys the power law 



/(ej ~ £ 



\l(q-\) 

k ' 



(22) 



as desired. On the other hand, for 0<^<1, /(e<.) has a cut-off at £,,^^ = 
[(l-^)jS*] in view of eq. (6). 

Since the number of degrees of freedom of the heat bath is assumed to be very large, 
it is appropriate to consider the relative probability of finding the system I as the ratio of 
the probability in the state with energy e^, relative to the fixed value of energy £, : 






(23) 



which is rewritten as 



A^k'^^i) 



' 'a,,[E-e,)^^ 



In 

V V 



^n[E-£) 



J) 



= e. 



^w(^_^^) KK(^-gJ)-ln,K(£-g.))]j- (24) 



Performing the expansion with respect to £^ and £, and keeping the leading order 
terms, we obtain 



A^k\£)^%[-^\£k-£))- 



(25) 
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Here, it should be noted that, to derive the relative probability of this form, we have to 
resort to the first principle calculation shown in eq. (24) and not to simply substitute the 
expression given in eq. (20) into the first part in eq. (23). This is a manifestation of the 
nonextensive (nonadditive) structure of the measure used here. 

In conclusion, we have shown that macroscopic thermodynamics of equilibrium 
concerning two systems in contact can be realized in ways other than the well-known 
Gibbsian. This has been achieved by relaxing the assumption of extensivity of the 
measure in phase space in a completely consistent way. Thus, the connection between 
macroscopic thermodynamics and microscopic statistical thermodynamics is established 
for systems obeying power- law canonical distributions. It should also be noted that this 
demonstration is in conformity with Tsallis nonextensive statistical mechanics [11]. We 
stress that the particular choice of the pair of functions, ln^(.r) and e^ix), in the 
present work is just one among many other possibilities. Thus, nonuniqueness of 
canonical ensemble theory for macroscopic thermodynamics of equilibrium is revealed. 
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